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ABSTRACT 


Of  the  many  Free  Electron  Lasers  (FELs)  in  the  world  today,  most  are  big  and 
expensive.  This  is  true  across  the  wavelength  spectrum,  from  long  to  very  short.  In 
contrast,  the  EEL  facility  in  progress  at  NFS,  which  will  initially  operate  at  long 
wavelengths  and  at  electron  energies  of  only  a  few  MeV,  is  inexpensive  and  smaller. 
However,  longer  wavelengths  lead  to  more  diffraction,  which  may  result  in  beam 
spread  and  interaction  with  the  undulator  surfaces.  Anticipating  the  possibility,  in 
this  thesis  we  analyze  mathematically  the  free  space  Hermite-Gaussian  modes  of  the 
optical  beam,  and  then  compare  them  to  the  Hybrid  modes,  where  waveguide  plates 
control  diffraction  along  one  axis,  allowing  free  space  diffraction  along  the  other  axis. 
We  continue  the  analysis  of  the  relativistic  electron  beam,  co-propagating  with  the 
optical  wave  in  the  Hybrid  Mode,  to  dehne  new  operating  condition  for  the  FEL. 
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I. 


INTRODUCTION 


The  Free  Electron  Laser  (EEL)  uses  a  beam  of  relativistic  electrons  co-propaga- 
ting  with  an  optical  beam  in  a  periodic  magnetic  field.  This  magnetic  field  causes 
the  free  electrons  to  oscillate,  allowing  coherent  energy  exchange  between  the  electron 
beam  and  the  optical  beam. 

The  first  development  relevant  to  the  EEL  was  the  microwave  tube,  invented  in 
the  1930s,  which  utilized  a  beam  of  free  nonrelativistic  electrons  in  a  closed  microwave 
cavity.  Later,  atomic  and  molecular  lasers  that  employed  a  system  of  bound  electrons 
and  an  open  optical  resonator  were  devised  in  the  1960s.  By  using  a  beam  of  free 
relativistic  electrons  and  the  open  optical  resonator,  J.  M.  J.  Madey,  the  inventor  of 
the  EEL,  combined  the  attributes  of  microwave  tubes  of  the  1930s  and  the  atomic 
lasers  of  the  1960s  [1]. 

Today  there  are  many  EELs  to  be  found  in  the  world,  and  most  of  these  pro¬ 
duce  high  quality  laser  light  at  IR  wavelengths  for  scientific  research.  EELs  generating 
short  wavelengths,  UV  and  X-ray,  have  been  studied  extensively  during  the  past  four 
decades;  these  require  large  and  expensive  accelerators  with  electron  energy  on  the 
order  of  GeV.  In  this  thesis,  I  will  focus  on  long  wavelength  EELs,  which  require  only 
a  few  MeV  for  the  electron  beam  energy  and  much  simpler  technology  with  lower 
costs.  This  thesis  work  will  provide  a  basic  theory  for  designing  EELs  to  operate  at 
long  wavelengths  in  the  laboratory.  Since  NFS  is  starting  a  EEL  facility  that  will 
first  reach  lower  energies  of  only  a  few  MeV,  the  EEL  considered  here  could  be  the 
first  to  operate  at  NFS. 

The  thesis  is  organized  as  follows:  Chapter  II  describes  the  EEL  system,  briefly 
discussing  how  EELs  work  and  their  significant  attributes.  Chapter  III  describes  the 
basic  theory  of  EEL  operation  dealing  with  the  electron  behavior  and  interaction 
between  the  electron  beam  and  the  optical  field.  This  chapter  will  be  limited  to 
plane  waves  to  characterize  the  optical  field.  However,  EELs  actually  operate  in 
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Gaussian  wave  modes.  So  Chapter  IV  offers  the  solution  of  the  wave  equation,  with 
an  emphasis  on  Hermite  Gaussian  resonator  modes.  Then  Chapter  V  explores  FEL 
theory,  including  the  features  of  Hermite  Gaussian  resonator  modes.  Chapter  VI 
deals  with  Hybrid  (HG-Wavegnide)  resonator  modes.  These  are  important  to  analyze, 
because  in  order  to  operate  a  FEL  at  long  wavelengths,  we  may  need  to  use  waveguide 
structure  capable  of  confining  the  radiation  that  otherwise  spreads  due  to  diffraction. 
New  research  is  then  presented  in  Chapter  VH  to  snggest  an  FEL  theory  applied  to 
Hybrid  (HG-Wavegnide)  resonator  modes,  which  is  the  ultimate  goal  of  this  thesis. 
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II.  FREE  ELECTRON  LASER  SYSTEM 
DESCRIPTION 


The  Free  Electron  Laser  system  consists  of  several  different  elements  involving 
interesting  physics.  Figure  1  shows  the  major  components  of  an  FEL  system  such  as 
the  Jefferson  Lab  14kW  IR  FEL,  and  the  proposed  NFS  FEL.  There  are  two  basic 
configurations  of  an  EEL,  the  oscillator  and  the  amplifier,  which  will  be  examined 
in  the  following  section.  Most  of  the  components  shown  in  Figure  1  are  common  to 
both  designs. 

A  pulsed  electron  beam  is  produced  when  an  injector  sends  the  beam  into 
an  accelerator,  such  as  superconducting  linear  accelerator.  The  beam  of  electrons 
is  accelerated  to  near  the  speed  of  light  utilizing  radio- frequency  (RF)  fields.  The 
relativistic  electrons  enter  the  undulator  and  then  get  “wiggled”  to  produce  laser  light, 
bringing  about  coherent  emission.  In  this  process,  only  a  small  fraction  (typically  ~ 
1-2%)  of  the  electron  beam  energy  is  converted  to  optical  energy.  In  some  cases,  the 
electrons  are  then  re-circulated  through  the  accelerator  in  order  to  recover  most  of 
their  remaining  energy.  Finally,  the  electron  beam  enters  the  beam  dump. 


» 


RF  Power 

Liquid  Helium  Refrigerator 


Beam 

Dump 


Injector 


FEL  Oscillator 


Undulators 


FEL  Amolifier 


Superconducting  Electron  Accelerator 


Seed  Laser 


Figure  1.  FEL  main  system  schematic.  Note  the  two  configurations  identified  as 
the  “Oscillator”  or  the  “Amplifier.”  Red  dots  represent  the  electron  pulses,  and  blue 
ellipses  describe  the  optical  pulses.  Erom  [2]. 


3 


A.  CONFIGURATIONS 


In  an  oscillator  configuration,  the  radiation  is  stored  and  amplified  between  two 
refiective  mirrors  as  shown  in  Figure  2.  The  optical  beam  is  then  outcoupled  through 
a  partially  transmitting  mirror.  Coherence  is  developed  through  many  passes  of  the 
optical  field  through  the  undulator. 


Figure  2.  FEL  schematic  in  an  oscillator  configuration.  From  [3]. 


In  an  amplifier  configuration,  the  optical  resonator  does  not  exist;  instead,  the 
system  employs  a  “seed  laser,”  which  is  an  external  laser  source.  Since  only  one  pass 
is  used,  the  amplifier  uses  a  longer  undulator  to  attain  more  gain.  Coherence  in  this 
case  is  established  by  the  seed  laser.  Figure  3  shows  a  typical  amplifier  configuration. 


Hlectfoa  Beam 


Figure  3.  FEL  schematic  in  an  amplifier  configuration.  From  [3]. 


B.  ATTRIBUTES 

The  FEL  has  some  unique  attributes  such  as  tunability,  designability,  high 
power,  efficiency  and  reliability.  Tunability  separates  the  FEL  from  a  conventional 
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laser  in  which  the  wavelength  is  determined  by  atomic  energy  levels,  n  an  FEL,  the 
resonance  condition  between  the  Lorentz  factor  7  of  the  electron  beam,  the  undu- 
lator  period  Aq,  and  the  dimensionless  undulator  parameter  K  determine  the  laser 
wavelength  A: 

^  _  Ao  (1  +  K^) 

V  ’ 

where  K  is  proportional  to  the  magnetic  field  of  the  undulator.  This  relationship 
means  that  by  changing  the  spacing  of  the  magnets,  the  strength  of  the  magnetic 
field,  or  the  energy  of  the  electron  beam,  the  FEL  can  be  continuously  tuned  to  dif¬ 
ferent  wavelengths.  Also,  FELs  have  been  designed  over  a  wide  range  of  wavelengths, 
from  microwaves  to  X-rays;  that  is,  they  possess  “designability.”  The  absence  of  a 
vulnerable  medium  such  as  gasses,  fluid  chemicals  or  solid  materials  means  that  the 
FEL  has  the  potential  to  achieve  extremely  high  power.  As  a  result,  the  FEL  can 
have  good  efficiency  (~  10%)  when  recirculating  the  electron  beam.  The  final  advan¬ 
tage  of  the  FEL  is  reliability:  existing  FEL  systems  now  run  twenty-four  hours  a  day 
for  weeks.  On  the  down  side,  the  FEL  is  as  huge  as  a  laboratory  and  considerably 
expensive. 

C.  COMPONENTS 

The  major  components  of  an  FEL  system  are  an  injector,  an  accelerator,  an 
undulator,  a  resonator  for  the  oscillator,  a  seed  laser  for  the  amplifier,  and  a  beam 
dump.  There  are  many  other  constituents  that  allow  the  FEL  to  operate,  such  as 
bending  magnets  for  steering  the  path  of  the  electron  beam,  an  optical  beam  transport 
system,  radiation  shielding  and  vibration  control. 

1.  Injector 

The  injector  generates  and  accelerates  the  beam  of  free  electrons  to  about  ~  5 
MeV  energy.  There  are  several  ways  to  produce  the  beam  of  free  electrons.  Typically, 
the  free  electrons  are  emitted  into  vacuum  through  either  thermionic  emission  or  the 
photoelectric  effect.  In  the  thermionic  cathode,  the  electrons  are  emitted  from  the 
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surface  of  the  cathode  when  the  cathode  is  heated  sufficiently.  A  photocathode  uses 
an  incident  drive  laser  to  excite  the  electrons  in  a  cathode.  By  using  the  photoelectric 
effect,  electrons  are  ejected  as  long  as  the  photon  energy  is  greater  than  the  cathode 
work  function.  The  released  electrons  can  be  emitted  into  a  Radio  Frequency  (RF) 
cavity,  that  accelerates  the  electrons  to  the  desired  energies  5  MeV)  before  they 
enter  the  main  accelerator.  The  electron  beams  usually  come  out  in  pulses  with 
excellent  quality  (a  narrow,  collimated  beam  with  low  energy  spread);  not  maintaining 
such  beam  quality  can  affect  FEL  system  performance.  An  RF  injector  is  shown  in 
Figure  4. 


Figure  4.  Schematic  of  an  injector.  From  [4]. 

2.  Accelerator 

After  the  injector,  the  electron  beam  is  now  passed  on  to  the  accelerator. 
The  accelerator  increases  the  energy  of  the  electron  beam  from  ~  5  MeV  to  the 
desired  level.  Direct  Current  (DC)  acceleration  is  possible,  but  most  FELs  use  a 
Radio  Erequency  Linear  Accelerator  (RE  LINAC)  to  obtain  a  final  energy  of  ~  100 
MeV.  The  accelerator  may  consist  of  several  evacuated  superconducting  metal  RF 
cavities.  The  electromagnetic  field  generated  in  each  cavity  should  be  in  phase  with 
the  injected  electron  pulses  to  provide  net  positive  acceleration. 
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3. 


Undulator 


The  undulator  is  a  significant  part  of  the  PEL  located  inside  the  optical  res¬ 
onator  cavity  and  typically  consists  of  alternating  permanent  magnets.  When  the 
electrons  leave  the  accelerator,  they  enter  the  periodic  magnetic  field  and  are  made 
to  wiggle  in  a  transverse  direction  by  the  Lorentz  force.  The  wiggling  motion  causes 
the  electrons  to  radiate  photons  and  amplify  the  fight  already  in  the  undulator.  As 
the  electrons  and  light  interact  in  the  undulator,  there  is  an  exchange  of  energy:  some 
electrons  gain  energy  from  the  fight  and  some  lose  energy  to  the  fight.  The  electrons 
that  gain  energy  move  a  little  bit  faster  along  the  axis  of  the  undulator  and  those  that 
lose  energy  move  a  little  bit  slower;  this  causes  “micro  bunching”  of  the  electrons  on 
the  scale  of  an  optical  wavelength.  Undulators  typically  have  either  linear  or  helical 
polarization.  A  linear  undulator  drives  the  electrons  to  move  sinusoidally  in  a  plane 
while,  a  helical  undulator  drives  them  to  a  corkscrew  motion.  A  schematic  of  an 
undulator  and  the  electron  beam  path  is  shown  in  Figure  5. 


Figure  5.  Schematic  of  an  undulator.  The  red  line  is  the  electron  beam  wiggling  along 
the  transverse  axis.  The  yellow  fines  indicate  the  stimulated  emission  propagating 
along  an  axis  of  forward  direction.  The  undulator  consists  of  magnets  alternating 
north  and  south  poles,  resulting  in  the  periodic  magnetic  fields.  From  [5]. 
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4. 


Resonator 


In  the  oscillator  configuration,  the  resonator  is  the  component  in  which  the 
light  is  stored  and  amplified  over  many  passes.  It  is  made  up  of  an  evacuated  cavity 
bounded  by  two  mirrors.  One  of  the  mirrors  has  high  reflectivity,  while  the  other 
allows  a  percentage  to  come  out  through  a  partially-transmitting  mirror.  The  light 
derived  from  this  out-coupling  mirror  is  the  usable  laser  light  of  the  FEL  system.  It  is 
obvious  that  for  maximum  coupling  efficiency,  electron  pulses  must  be  synchronized 
with  the  bounce  time  of  optical  pulses  between  resonator  mirrors.  A  schematic  of  the 
electron  and  optical  pulses  passing  through  the  undulator  in  the  resonator  is  shown 
in  Figure  6. 


dldctron  beam  pulsas 

wiggler 

f  •  •  • 


100  % 


reflactlng 

mirror 


• 

partially 

transmitting 

mirror 


optical 

laser 

pulses 


Figure  6.  Schematic  of  a  resonator.  From  [2]. 


5.  Beam  Dump 

The  beam  dump  is  the  last  part  of  the  electron  beam  path.  It  consists  of  an 
absorbing  material  (such  as  a  large  metal  block)  surrounded  by  radiation  shielding. 
Extracting  energy  from  the  electron  beam  amplifies  the  optical  beam,  but  the  extrac¬ 
tion  of  energy  is  typically  only  a  few  percent,  leaving  a  high-energy  electron  beam 
which  results  in  poor  system  efficiency,  and  lots  of  potentially  dangerous  radiation 
emitted  at  the  dump.  If  the  system  uses  recirculation  of  the  electron  beam,  the  elec¬ 
trons  enter  to  the  LIN  AC  a  second  time  180°  out  of  phase,  decelerating  electrons, 
and  giving  most  of  their  energy  back  to  the  RF  field.  The  recirculated  electron  beam 
has  approximately  ^3  MeV  average  energy  at  the  end  of  the  LINAC,  having  lost  a 


few  percent  to  the  laser  beam  and  leaving  about  the  same  energy  that  the  electron 
beam  had  when  it  was  generated  at  the  injector.  This  greatly  improves  overall  system 
efficiency,  and  reduces  radiation  shielding  requirements  at  the  dump.  A  schematic  of 
a  beam  dump  is  shown  in  Figure  7. 
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Figure  7.  Schematic  of  a  beam  dump.  From  [7]. 
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III. 


BASIC  FREE  ELECTRON  LASER  THEORY 


We  will  now  explore  the  mathematics  that  describe  the  operation  of  the  PEL. 
The  physics  of  the  electron’s  behavior  and  the  interaction  between  the  electron  beam 
and  optical  held  are  pivotal  to  the  operation  of  EEL  systems.  In  the  following  sec¬ 
tions,  the  X  and  y-sods  designate  the  transverse  axes  while  the  .2- axis  designates  the 
longitudinal  axis  passing  through  the  center  of  the  undulator. 

A.  ELECTRON  BEHAVIOR 

The  electrons  in  the  undulator  have  a  microscopic  evolution. 

1.  Undulator  Field  and  Transverse  Motion 

Exploring  the  behavior  of  electrons  in  the  magnetic  held  of  the  undulator  is 
necessary  to  understand  the  physics  of  the  EEL  interaction.  We  will  deal  with  rela¬ 
tivistic  electrons  entering  a  helical  undulator  held  (which  is  simpler  mathematically 
than  a  linear  undulator)  in  helical  orbits  around  the  ^-axis  in  the  presence  of  the 
optical  held.  An  electron  passing  through  the  undulator  is  exposed  to  three  helds: 
the  magnetic  held  from  the  undulator  magnets,  the  magnetic  held  of  the  optical  wave 
and  the  electric  held  of  the  optical  wave,  as  shown  in  Figure  8. 


Figure  8.  For  convenience,  undulator  and  Optical  Fields  in  the  case  of  linear  undula¬ 
tor.  From  [2] 
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The  magnetic  field  of  the  helical  undulator  is 


Bu  =  B  [cos  {koz) ,  sin  (koz) ,  0] 


(iii.i) 


where  B  is  the  magnetic  field  amplitude,  fco  =  27r/Ao  is  the  undulator  wave  number 
and  Ao  is  the  undulator ’s  period  which  is  the  distance  of  one  complete  cycle  of  the 
magnetic  field  in  the  undulator.  The  corresponding  electric  and  magnetic  field  of  the 
optical  field  are 

Er  =  £"  (cosi/i,  —  sin-^,  0)  (IIL2) 

Br  =  E  (sin  'll;,  cos  'ip,  0)  (III. 3) 


where  E  is  the  optical  wave  amplitude  in  cgs  units,  'ip  =  kz  —  wt  +  (p,  k  =  27: / \  is  the 
optical  wave  number,  w  =  kc  is  the  optical  angular  frequency,  and  cp  is  the  optical 
phase. 

The  relativistic  Lorentz  force  equations  determine  the  electron  motion  in  the 
undulator.  The  relativistic  Lorentz  force  equations  in  cgs  units  [1]  for  an  electron  are 


4  (7^)  =  -—  [Er  +  /3  X  (B„  +  Br)]  , 
dr  me 


d7 

dt 


e 

me 


/3-Er, 


7  ^  =  1  - 


(111.4) 

(111.5) 

(111.6) 


where  v  =  /3c  is  the  electron  velocity,  e  =|  e  |  is  the  electron  change  magnitude,  m  is 
the  electron  mass,  c  is  the  speed  of  the  light,  and  7  is  the  relativistic  Lorentz  factor. 
Substituting  the  fields  from  Equations  (III.I),  (HI. 2),  and  (IIL3)  into  the  (IIL4),  we 
get 

{ll3x)  =  — —[Ecos'ip{l-  /3^)  -  f3^B sin  (koz)] ,  (111.7) 

dr  me 

=  — —  [-B  sin  V-  (1  -  (3z)  +  l3zB  cos  (koz)] ,  (III. 8) 

cir  TTT/e 

d  6 

i'yf^z)  = - [E  {(ix  cos  pJ  -  Py  sin  p)  +  B{Px  sin  {k^z)  -  Py  cos  (koz)}] .  (III. 9) 

dr  TThe 
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Since  the  electrons  are  traveling  close  to  the  speed  of  light  along  the  2;  axis,  ~  1, 
so  that  terms  containing  (1  —  l3z)  can  be  ignored  in  the  transverse  x-y  components  of 
Equations  (III. 7)  and  (IIL8).  In  other  words,  E  {1  —  Pz)  is  small  compared  to  PzB, 
which  means  magnetic  fields  derived  from  an  undulator  have  a  bigger  effect  on  the 
electron’s  transverse  motion  than  the  optical  fields.  Integrating  Equations  (HI. 7)  and 
(III. 8)  with  respect  to  time  and  setting  the  constant  of  integration  to  zero,  indicating 
perfect  injection  into  helical  orbits,  we  get 


—  cos  {hz)  , 

7 

(III.IO) 

—  sin  {k^z) , 

7 

(III.ll) 

where  K  =  ei?rmsAo/27rmc^  is  the  dimensionless  undulator  parameter  and  Brms  is 
the  rms  value  of  the  undulator  field.  Eor  the  helical  undulator  Brms  =  B;  for  the 
linearly  polarized  undulator  Brms  ~  Bj'/^.  Integrating  once  again  and  using  k^z  — 
koVzt  =  koPzCt  ~  koct  =  ujot,  we  get 

rK  K\n 

X  (t)  = - sin  (uot)  = - sin  (uiot) ,  (III.  12) 

(UqT  27r7 

cK  KXn 

y  (t)  = - cos  (uot)  — - cos  (uot) ,  (III.  13) 

27r7 

which  describes  the  transverse  (helical)  motion  of  the  electron. 


2.  Longitudinal  Motion  and  the  Electron  Phase 

Now  consider  the  effects  of  the  optical  electric  field.  We  scrutinize  the  micro¬ 
scopic  electron  motion  along  the  longitudinal  axis  which,  while  small  in  magnitude, 
is  key  to  EEL  operation.  Substituting  Equations  (III. 10)  and  (III. 11)  into  (III. 9)  and 
using  a  trigonometric  identity  gives 

Ipz  =  cos  {koZ  +  p)-  ppz-  (III.  14) 

'jmc 

In  order  to  get  7,  substitute  Equations  (III. 10),  (III. 11)  and  (III. 2)  into  (III. 5)  and 
get 


dq 

dt 


=  7 


eKE 

jmc 


cos  {koz  +  p)  ■ 
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(III.15) 


Substituting  Equation  (111.15)  into  (111.14),  we  arrive  at 


•  r  ±\  H/ 

=  — cos  (fco^  +  V')  •  (III. 16) 

'y^mc 

Also,  we  can  rewrite  the  Equation  (III.  15)  by  using  the  definition  of  ip  and  the  “elec¬ 
tron  phase,”  (p  =  {k  +  ko)  z  —  ut,  which  describes  the  electron’s  phase  with  respect  to 
the  combined  optical  and  undulator  fields,  to  obtain 

7  =  cos  ((  +  (p) .  (III. 17) 

qmc 

Since  ko,  k  and  u  are  fixed,  the  dynamical  variable  in  ^  (t)  is  .2  (t)  which  represents 
the  electron’s  microscopic  position,  (p  oc  kz  (t).  Note  that  when  cos  {(p  +  cp)  >  0,  the 
electron  energy  increases  by  absorption,  7  >  0;  when  cos  {C  +  <P)  <  0,  the  electron 
energy  decreases  by  stimulated  emission,  7  <  0.  This  causes  the  electron  beam  to 
bunch  and  to  radiate  coherently. 

3.  Resonance  Condition  and  Pendulum  Equation 

EEL  resonance  occurs  when  the  energy  exchange  between  the  electron  and  the 
optical  beam  is  optimum.  The  electron  travels  down  the  z-axis  with  a  speed  while 
the  optical  pulse  travels  at  the  speed  of  light,  c.  The  relationship  between  these  two 
speeds  brings  about  the  “electron-photon  race,”  which  is  illustrated  in  Figure  9. 


Figure  9.  Electron-photon  race.  From  [2]. 

The  time  for  an  electron  to  travel  through  one  undulator  period  is  Xq/(3zC-  During  that 
time,  one  wavelength  of  light  is  emitted  because  the  electron  executes  one  oscillation 
in  the  undulator.  The  wavelength  of  light  then  travels  a  distance  A  -|-  Aq  at  speed  c. 
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So, 


-^0 


A  +  Aq 


From  (111.18)  we  see  that  the  laser  wavelength  at  resonance  is 

Ao  (1  —  Pz) 


A  = 


Pz 


(III.18) 


(III.19) 


Now  we  hnd  an  expression  for  Pz  in  terms  of  the  dimensionless  undulator 
parameter,  K.  From  Equations  (III.  10)  and  (III.  11),  we  find  that  P±^  = 
Therefore,  the  relativistic  Lorentz  factor,  7,  which  indicates  the  relationship  between 
the  electron  energy  and  the  z-axis  velocity  cPz  can  be  expressed  as 


=  1  -  /3^  =  1  -  =  1  -  —  -  P/. 

y2 

For  relativistic  electrons,  7  1,  we  can  use  the  binomial  expansion  to  get 

1  +  /F2 


(III.20) 


Pz 


1  - 


272 


(III.21) 


We  can  rewrite  the  resonance  condition,  Equation  (III. 9),  in  terms  of  K  to  get 

Ao  (1  +  K^) 


A 


272 


(III.22) 


By  changing  the  physical  parameters,  such  as  Aq,  K  and  7,  we  can  get  a  wide  range  of 
wavelengths  of  the  optical  beam,  which  indicates  the  EEL’s  tunability.  Also  we  can 
now  see  that  the  wavelength  of  the  optical  beam  A  is  much  smaller  than  the  undulator 
period  Aq  because  of  the  small  factor  1/7^. 

Now  we  are  ready  to  derive  the  EEL  “pendulum  equation”  by  combining  some 

o 

of  the  previous  equations.  Let  us  introduce  an  operator  (..)=  d(..)/  dr  which  can 
be  expressed  as 


_d_ 

dr 


d  dt 


L_d_ 
c  dt  ’ 


dt  dr 

where  r  =  ct/L  is  the  “dimensionless  time,”  where  L  is  the  undulator  length,  so  that 
r  =  0  — >  1  along  the  undulator.  By  using  this  relation,  we  can  determine  the  rate  of 
change  of  the  electron  phase  along  the  undulator  as  the  electron  phase  velocity, 

dC 


u  —  — C—  L  [{k  +  ko)  Pz  —  k] , 


(111.23) 


15 


where  the  derivative  of  C  (electron  phase)  is  with  respect  to  r;  again  L  =  NXq  is  the 
length  of  the  undulator,  so  that  0  <  r  <  1.  By  substituting  Equation  (III. 21)  into 
(III. 23)  and  using  k  ^  ko  near  resonance,  we  arrive  at  [6] 

dn  =  AttN—,  (III.24) 

7 


relating  changes  in  beam  energy  to  the  electron  phase  velocity.  Combining  Equation 
(III.  15)  and  (III. 24)  gives 


dn  L  du  L  diy  dq 
dr  c  dt  c  dq  dt 


LdTTNeKE 

- cos  (C  +  0) , 

c  q  qmc 


C=  I  a  I  cos(C  +  ^) , 


(III.25) 


where  |  a  \=  dirNeKLE/^'^mc^  is  the  “dimensionless  optical  field  amplitude.”  Equa¬ 
tion  (III. 25)  has  the  form  of  the  pendulum  equation,  and  describes  the  microscopic 
motion  of  the  electrons  in  phase  space,  (^,  v). 


4.  Phase  Space  Evolution 

Since  phase  velocity  is  associated  with  the  electron  energy  according  to  Equa¬ 
tion  (III. 24),  phase  space  diagrams  show  an  energy  exchange  between  electrons  and 
the  optical  field.  Assuming  |  a  \  and  (p  are  approximately  constant,  the  electrons 
follow  phase  space  paths  given  by  [8] 


+  2  \a  \  [sin  (C  +  0)  -  sin  (Co  +  (p)] ,  (III. 26) 


where  Co  and  Pq  are  the  initial  electron  coordinates  in  phase  space.  A  special  phase 
space  path  is  the  EEL  separatrix  [9],  starting  at  the  unstable  fixed  point  (Co  =  37r/2 
and  pq  —  0),  so  that 

p'^  =  2  \  a  \  [1  d-  sin  (C  +  (p)]  ■  (III. 27) 

Erom  this  equation,  we  can  see  that  the  peak  to  peak  height  of  the  separatrix  is 
4|  a  1^^^.  A  larger  |  a  \  results  in  the  phase  space  paths  evolving  much  faster  and  a 
separatrix  increased  in  height.  Inside  the  separatrix,  the  electrons  follow  closed  orbits 
(corresponding  to  a  pendulum  swinging  back  and  forth).  Outside  the  separatrix,  the 
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electrons  follow  open  orbits  (like  a  pendnlnm  swinging  over  the  top).  Figure  10  shows 
the  evolution  of  20  sample  electrons  at  resonance  t'o  =  0  and  weak  field  ao  <  tt  during 
the  dimensionless  time  from  r  =  0  (yellow  dots)  to  r  =  l(red  dots). 

0.2 

0 

0.1 

0 

-1t/2  ^  31C/2  0  t  1 


Figure  10.  Phase  space  evolution  of  electrons  injected  at  resonance  {uq  =  0),  and 
corresponding  optical  field  gain  and  phase  evolution.  From  [4]. 


Figure  11.  Phase  space  evolution  of  electrons  injected  off  resonance  {pq  =  2.6),  and 
corresponding  optical  field  gain  and  phase  evolution.  Periodic  boundary  conditions 
are  applied  to  keep  the  electron  phase  (  between  — 7r/2  and  37r/2.  From  [4]. 


The  symmetrical  distributions  of  the  electrons  indicates  that  about  half  of  the  elec¬ 
trons  gain  energy  from  the  optical  field  (paths  shift  up)  and,  the  other  half  lose  energy 
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to  the  optical  field  (paths  shift  down) .  We  can  see  a  bunching  of  electrons  occurs  but 
there  is  no  gain  since  an  equal  number  of  electrons  gained  and  lost  energy. 

If  electrons  were  injected  slightly  off  resonance  {uq  =  2.6),  as  shown  in  Figure 
11,  more  electrons  would  shift  to  a  lower  value  of  phase  velocity  than  shift  to  higher 
energy.  There  is  an  obvious  bunching  of  electrons  near  ^  =  tt,  which  causes  the 
positive  gain. 

Gain  is  the  ratio  of  the  change  to  the  initial  value  in  the  power  of  the  optical 
field  as  it  passes  along  the  undulator.  We  can  express  the  gain  as 

a^(T)  —  ao^ 


G{t)  = 


ao" 


(III.28) 


where  a(r)  is  the  field  along  the  undulator  at  r  and  ao  is  the  initial  field  starting  at 
beginning  of  the  undulator. 


B.  ELECTRON-OPTICAL  FIELD  INTERACTION 

We  derived  the  pendulum  equation  which  represents  the  microscopic  evolution 
of  the  electrons  in  the  FEL  interaction  region.  We  now  need  to  derive  an  equation 
describing  the  self-consistently  evolving  optical  field. 

1.  Maxwell’s  Equations 

Now,  consider  the  interaction  between  the  electron  beam  and  optical  field,  and 
start  with  Maxwell’s  equations  in  cgs  units  [10]: 

V  •  E  =  dyrp 


VB  =  0 

V  X  E  =  —  — 
c  ot 


dyr  1 5B 

V  X  B  =  — J  +  -  — , 
c  c  ot 


(III.29) 

where  p  is  the  charge  density,  J  is  the  current  density.  From  Maxwell’s  equations,  we 
can  derive  the  full  wave  equation  associated  with  transverse  current  [8], 

^'-^^)A(xG)  =  --Jx(x,t),  (III.30) 
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where  A(x,  t)  is  the  optical  vector  potential,  c  is  the  speed  of  light,  and  J_l(x,  t)  is  the 
transverse  current  density;  in  our  case,  due  to  the  wiggling  motion  of  the  electrons 
traveling  through  the  undulator. 

2.  FEL  Wave  Equation 

For  our  helical  undulator,  the  optical  vector  potential  can  be  written  as 

A(x,t)  =  (III.31) 

k 


where  E  =\  E  \  is  the  complex  laser  electric  field,  a  =  kz  —  ujt  is  the  carrier 

wave,  i  =  (— ^1,0)  is  the  laser  field’s  polarization  vector  for  the  case  of  circular 
polarization  from  a  helical  undulator,  and  x  =  xi  +  r/j  +  zk  is  the  position  vector. 
By  substituting  Equation  (III.  31)  into  (III.  30)  and  assuming  that  the  optical  wave 
amplitude  |  E{'x,t)  \  and  phase  0(x,  t)  are  slowly  varying  along  the  2;  axis,  we  arrive 
at  [2] 


V I  +  2ik  (  —  I  — 
oz  c  ot 


E  =  Jx  -  re-*". 


(III.32) 


which  is  the  parabolic  or  paraxial  wave  equation  with  a  source  current  J^.  When  we 
change  the  variables  to  a  reference  frame  that  “follows  the  light,”  u  =  z  —  Lt  =  z  —  ct 
and  T  =  ctjL,  then  we  can  rewrite  the  operator  [11] 


d  Id  Id 
dz~^  cdt  L  dr 


(III.33) 


The  source  current  which  is  the  sum  of  all  single-particle  currents  can  be  expressed 
as 

Jj_  =  —ec  _Lh^^^(x  -  Xj(t)),  (III. 34) 

i 

where  (5*^^)(x  —  Xj(t))  =  S{x)S{y)S{z)  is  the  three-dimensional  Dirac  delta  function, 
and  Xj(t)  is  the  electron  position  at  time  t.  The  transverse  motion  of  an  electron 
l3j_  which  contributes  to  the  transverse  current  in  the  undulator  can  be  expressed 
as 

I3j_  =  — —  (cos  A:o.2,  sin  A:o.2,  0)  =  Re  (  .  (III. 35) 

7  \  7  / 
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Substituting  Equation  (III.  35)  into  (III.  34),  and  (III.  34)  into  (III.  32)  gives 


2ik  d 

~TdT 


E{v,t)  =  —A.meKkp{-x.,t)  {e  , 


(III.36) 


where  r  =  xi+yj,  and  p(x,  t)  =  f  (x  —  Xj(t))  dV  is  the  local  electron  density  in 

a  small  volume  element  dV,  and  <  ...  >  indicates  an  average  over  sample  electrons  in 
the  volume  element.  If  we  multiply  Equation  (III.  36)  by  {A7:NeKL/j‘^mc^){—Li/2k) 
and  define  the  complex  laser  field  a  =  |  o  |  with  the  previous  laser  field  amplitude 


a  1=  4:7rNeK LE / ,  we  can  get 


2  d 

2k  ^  dr 


a 


(r,T)  =  -  (je  , 


(III.37) 


where  j  =  p/^^m(?  is  the  “dimensionless  current  density”  which  repre¬ 

sents  another  important  EEL  parameter  (maybe  the  most  important).  By  introduc¬ 
ing  the  dimensionless  transverse  coordinates  x  =  x{k/2L'f^^ ^  y  =  y{k/2L)^^‘^  and 
replacing  the  Laplacian  operator  with  a  dimensionless  Laplacian  operator,  V_l^  = 
{k/2L)'V']_,  we  arrive  at  the  EEL  wave  equation  [11] 

d 


rVl  + 


dr 


i{r,T)  =  -  {je  , 


(III.38) 


where  +  dy^,  and  the  coordinates  {x,y,z,T)  are  all  now  dimensionless. 

When  diffraction  is  neglected,  the  EEL  wave  equation  can  be  written  in  the  simplest 
form 

a=  -j  (e-*^)  .  (III.39) 


The  dimensionless  current  density  j  measures  the  coupling  between  the  electron  beam 
and  the  optical  mode,  and  measures  the  amount  of  electron  bunching.  When  j 

is  large  {j  S>  tt),  the  coupling  is  large  and  we  have  high  gain;  when  j  is  small(j  <  tt) 
we  have  low  gain. 


3.  Gain 

Until  now,  we  have  stndied  the  interaction  between  free  electrons  in  the  beam 
and  laser  light  through  the  Lorentz  force  equations  in  the  helical  undulator.  Previ- 
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ously,  we  have  understood  the  bunching  of  electrons  in  phase  space  resulting  in  coher¬ 
ent  radiation  through  the  pendulum  equation.  Now,  we  have  coupled  this  bunching 
to  the  laser  beam  through  the  PEL  wave  equation. 

When  we  assume  a  weak  field  |  o  tt  and  low  gain  j  <  tt,  we  can  apply  a 
pertubation  expansion  to  the  pendulum  equation,  resulting  in  PEL  gain,  [2] 

2  —  2cos(r^o'^)  “ 

where  a  filling  factor  F  represents  a  ratio  of  cross  sectional  areas  of  electron  beam  to 
optical  mode  such  as  F  =  {rh/woY  at  the  waist  of  the  mode,  where  is  the  electron 
beam  radius  and  wq  is  the  optical  mode  waist  radius.  Gain  is  proportional  to  the 
dimensionless  current  density  j  and  filling  factor  F,  and  the  final  gain  spectrnm  at 
r  =  1  depends  only  on  the  initial  electron  phase  velocity  Pq  as  shown  in  Eigure  12. 
As  seen  before,  there  is  no  gain  at  resonance,  pq  =  0,  and  the  peak  of  G/jF  (gain 
normalized  to  j  and  F)  is  abont  13.5%  and  occnrs  at  pq  ~  2.6.  / / 


(III.40) 


G(r)  = 


q('^) 

Oo^ 


-1=  jF 


Eigure  12.  Single  pass  G/jF  versus  initial  electron  phase  velocity  pq  at  t  =  1. 
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4. 


Extraction 


In  order  to  describe  the  energy  exchange  between  the  electron  beam  and  optical 
field,  we  estimate  the  extraction  r]  defined  as  the  fractional  energy  transferred  from 
the  electron  beam  to  the  optical  field  during  a  single  pass  through  the  undulator  [6] : 


V  = 


_  Eopt 


(AE)  (Ay) 


(III.41) 


Eb  Eb  7 

where  Eopt  and  Eb  are  the  energy  in  the  optical  field  and  the  electron  beam,  respec¬ 
tively,  in  a  common  volume  element  of  the  beams,  and  AE  =  Aymc^  is  the  net 
change  in  the  electron  beam  energy  through  the  undulator.  Applying  Equation  (III. 
24)  to  (III.  41)  allows  us  to  have  the  following  result: 


_  (Ay)  _  {Ap)  _  J_ 
y  dvrA  ~  2N 


(III.42) 


where  it  has  been  assumed  that  approximately  half  the  electrons  lose  energy  across 
the  positive  gain  bandwidth  Ap  ~  tt  illustrated  in  Figure  12.  We  know  the  height  of 
the  separatrix  is  Ay/a,  as  presented  below  Equation  (III.  27).  This  height  corresponds 
to  a  change  of  an  electron  phase  velocity  Ap,  when  an  electron  is  in  the  largest  closed 
orbit.  Suppose  that  only  half  of  the  electrons  are  in  closed  orbits,  then  (Ai^)  would  be 
2^/a.  Now,  we  need  to  prove  that  a  ~  by  giving  a  short  account  of  “Synchrotron” 
oscillation. 

Another  argument  starts  from  phase  space  considerations.  Start  with  the 
pendulum  equation  with  low  gain  j  <  vr  and  strong  fields  |  n  tt.  Assuming  |  a  |~ 
constant  and  0  =  0  yields  the  pendulum  equation  as 


c=  I  a  I  cos  (C) . 


(III.43) 


For  an  electron  trapped  in  strong  field  near  the  stable  fixed  point  (0o  =  7r/2,  0)  in  the 
phase  space,  its  motion  can  be  described  by 


C(^)  =  ^  + 


(III.44) 
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where  Co  =  '7i'/2  is  the  initial  electron  position  and  x{t)  -C  tt  [2],  Substituting  Equa¬ 
tion  (III.  44)  into  (III.  43)  gives 


oo  I 

x  =  \  a 


a  \  X  =  —p.  X, 


(III.45) 


where  Ug  is  the  synchrotron  frequency  about  the  stable  fixed  point  Co  =  7r/2.  We 
recognize  that  the  simple  harmonic  oscillation  of  electrons  around  the  stable  fixed 
point  (7r/2,0)  has  the  frequency  Ug  =  \  d  The  trapped  electrons  will  execute 
approximately  half  a  synchrotron  oscillation  at  saturation,  or  ~  tt  when  |  a  |~  tt^. 
Assuming  about  half  the  electrons  are  trapped,  we  arrive  at  the  estimate  for  extraction 


_  (At)  ^  (At/)  _  2^  _ 

7  4ttN  4ttN  2N 


(III.46) 
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IV. 


HERMITE  GAUSSIAN  RESONATOR  MODES 


In  the  previous  chapter,  we  examined  the  interaction  between  the  electron 
beam  and  the  optical  field  inside  the  undulator.  This  chapter  examines  the  wave 
nature  of  laser  beams  as  they  propagate  in  free  space.  Up  to  this  point,  we  have 
assumed  the  optical  beam  to  consist  of  plane  waves  neglecting  diffraction.  The  wave 
usually  takes  the  form  of  a  Gaussian  beam  as  defined  by  the  optical  cavity  bounded  by 
spherical  resonator  mirrors.  In  order  to  understand  the  characteristics  of  a  Gaussian 
beam,  we  will  find  solutions  to  the  free-space  wave  equation  including  diffraction. 


A.  SOLUTION  OF  WAVE  EQUATION 

The  characteristics  of  Gaussian  laser  beams  in  a  resonator  have  been  analyzed 
by  Siegman  [12],  Kogelnik  [13],  Pedrotti  [14],  and  others.  This  chapter  follows  their 
analysis. 


1.  Scalar  Wave  Equation 

Let  us  consider  a  general  wave  equation  for  the  vector  potential  A  without  a 
current  source, 


=  0  or  VA=i^.  (IV.l) 

Assuming  A{x,y,z,t)  =  u{x,y,  z)T(t)  and  utilizing  a  process  of  separation  of  vari¬ 
ables,  we  have 


1  11  d^T 

— =  —  —  — =  constant  =  —  k  >  0,  (IV. 2) 

A  ^  y  ^  dU  ^ 

space  time 

where  the  left-hand  side  is  only  a  function  of  space  and  the  right  side  is  only  a  function 
of  time,  so  that  u  represents  a  component  of  the  electromagnetic  field  or  of  the  vector 
potential  A.  When  we  consider  only  the  scalar  function  of  space,  u,  we  arrive  at  a 
scalar  wave  equation,  the  Helmholtz  equation  [15] 


(V^  -I-  k‘^)  u(x,  y,  z)  =  0, 


(IV.3) 
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where  k  =  27r/A  is  the  propagation  constant  in  the  medium,  with  a  harmonic  time 
variation, 

2.  Fundamental  Mode 

For  an  optical  beam  traveling  along  the  ^-direction.  Equation  (IV. 3)  becomes 

u{x,  y,  z)  =  V’(x,  y,  (IV.4) 


where  ^  is  a  slowly  varying  complex  valued  function,  and  describes  a  plane- 
wave  in  the  ^-direction.  Generally,  ip  represents  the  distinction  between  the  actual 
laser  beam  and  a  plane  wave,  which  we  will  examine  in  detail  below.  Substituting 
Equation  (IV.4)  into  the  scalar  wave  Equation  (IV. 3)  yields,  in  cartesian  coordinates, 
the  reduced  equation 


d'^ip  d'^ip  d'^ip  ,  dip 


(IV.5) 


As  ranges  over  an  optical  wavelength  the  change  in  ip  will  be  small,  also  the  variation 


of  Ip  with  respect  to  2;  will  be  much  smaller  than  with  respect  to  x  and  y.  This 
slowly  varying  dependence  of  on  2;  can  be  expressed  mathematically  by  the  paraxial 
approximation  [12]: 


< 

,  and 

dz'^ 

oz 

") 

dx"^ 

dy'^ 

(IV.6) 

We  can  then  discard  the  second  derivative  of  ip  with  respect  to  in  Equation  (IV.5) 
and  obtain  the  paraxial  wave  equation: 


d‘^ip  d‘^ip  dip 

_ L.  _i _ L.  _i_  Ooh——  =  n 

dx^  dy^^  dz 


(IV7) 


A  form  of  solution  for  Equation  (IV. 7)  in  the  fundamental  mode  can  be  written  as 
[13] 


ip{x,y,z)  =  exp 


P{z)  + 


k 


2q{z) 


(IV.8) 


where  =  x"^  +  y^,  is  a  complex  phase  shift  which  is  related  to  the  propagation 
of  the  optical  beam,  and  q{z)  is  the  complex  curvature  of  the  phase  front  which 
represents  the  spherical  shape  of  the  Gaussian  beam.  We  need  to  find  these  unknown 
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functions,  P{z)  and  q{z).  In  order  to  get  each  term  of  Eqnation  (IV. 7),  we  take 
derivatives  of  ip  with  respect  to  each  component.  We  then  have  the  following  sources 


d^ip 

dx^ 

d^ip 

w 


ik  9  ,  , 
- Ip, 

q  q‘^ 

ik  k"^  o  ,  , 

- -y^]  iP, 

q  q^ 


(IV.9) 
(IV.  10) 


and 


dip  _  ( .dP 
dz  \  dz 


ik 


Ip.  (IV.ll) 

Inserting  Eqnations  (IV.9),  (IV. 10),  and  (IV.ll)  into  (IV. 7)  and  dividing  throngh  by 


_  2^ 

2q‘^^  dz 


Ip,  we  arrive  at 


\dz  J  \dz  q  '  (  ^  ) 

In  order  for  Eqnation  (IV.  12)  to  be  valid  for  any  r^,  the  expressions  inside  both 
parentheses  must  be  zero.  Therefore,  we  obtain  the  relations 


^  =  1 

dz 


(IV.  13) 


and 


(IV.  14) 


dP  _  i 
dz  q 

At  this  point,  we  should  introduce  two  important  parameters  of  the  Gaussian 
beam,  w{z)  and  R{z),  where  w(z)  is  the  beam  radius  or  spot  size,  and  R{z)  is  the 
radins  of  curvatnre  of  the  wavefront  as  shown  in  Eignre  13.  These  parameters  are 
associated  with  the  waist  spot  size,  Wq,  and  the  ratio  zj Zq  by  the  formulas  [12] 

.2 


w{zy 


Wo 


1  +  ^ 

Zo^ 


(IV.15) 


which  can  be  derived  from  Equation  (IV.  13)  and 


R{z)  —  z  i  1  -\- 


Zo 


(IV.  16) 


where  zq  =  /\,  the  Rayleigh  length,  is  the  distance  which  the  beam  travels  from 

the  waist  to  where  the  beam  area  doubles. 
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Figure  13.  Notation  for  a  fundamental  Gaussian  beam  diverging  away  from  its  waist. 
Wq  is  the  waist  spot  size,  w{z)  is  the  beam  radius,  and  R{z)  is  the  radius  of  curvature. 
After  [12] 

We  now  consider  Equation  (IV.  13);  integrating  this  equation  gives 

q  =  z  +  qQ  =  z-izo,  (IV.17) 

where  the  initial  value  is  go  =  —mwQ^/\  [13].  Starting  with  Equation  (IV.17)  and 
using  Equations  (IV.  15)  and  (IV.  16),  we  get  the  following  result: 


1 

Q 


Z^  +  Zo^ 


+  i- 


^0 


.2^  +  .20^ 


1  A 

=  ^ - ; 


R 


TTW^ 


(IV.  18) 


(IV.  19) 


For  convenience  of  later  calculation,  we  can  rewrite  this  equation  as 

k  k  ikX  k  i 

2q  2R  2ttw‘^  2R  w‘^  ’ 

where  k  =  2tt/X.  To  determine  P{z)^  we  insert  Equation  (IV.17)  into  (IV. 14)  and  we 

have 

dP  iz  Zq 

dz  z^  +  z^  z^  +  z^ ' 


(IV.20) 
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Integration  of  this  equation  yields  the  following  result 


P{z)  =  i 


-1  ^ 

tan  — . 
zo 


By  using  Equation  (IV.  15)  we  can  rewrite  this  as 


(IV.21) 


iP(z)  =  In  f— )  -  —  ,  (IV.22) 

V  W  /  Zq 

where  the  tan“^  (z/zq)  represents  the  phase  shift  difference  between  the  Gaussian 
beam  and  a  plane  wave,  and  the  amplitude  factor,  wq/w,  describes  the  decrease  of 
the  beam  amplitude  on  the  axis  because  of  the  expansion  of  the  beam.  Substituting 
Equations  (IV.  19)  and  (IV.22)  into  the  trial  solution.  Equation  (IV. 8),  we  can  finally 
obtain  the  solution  of  the  fundamental  mode: 


w 


exp 


-1 


^0 


(IV.23) 


Using  this  solution,  Eqnation(IV.4)  can  be  written  in  the  form 

,  ^  ^ 

/  s  Wq  .^2/  2 

u[r,  z)  =  — e  '  exp 


Amplitude 


kz 


kr'^ 


—  tan  ^  — 
zo 


z 


2R 


V 


Phase 


/J 


(IV.24) 


Gaussian  beams  do  not  change  their  profile,  but  beams  which  have  long  wave¬ 
lengths  will  be  expanded  significantly  due  to  diffraction  as  they  propagate.  Eigure  14 
shows  well  the  expansion  of  Gaussian  beams.  In  order  to  show  the  expansion,  we  set 
the  wavelength  to  A  =  300pm. 


3.  Higher  Order  Modes 

In  the  previous  section,  we  discussed  only  the  lowest-order  solution  of  Equa¬ 
tion  (IV. 7).  There  are  other  higher-order  solutions  for  the  Equation  (IV. 7);  these 
can  take  the  form  either  of  Hermite-Gaussian  functions  in  cartesian  coordinates  or 
of  Laguerre-Gaussian  functions  in  cylindrical  coordinates  [12].  The  properties  of 
Laguerre-Gaussian  modes  in  the  EEL  were  analyzed  by  A.  Kampouridis  [4] .  Through¬ 
out  the  rest  of  this  thesis,  we  will  use  the  cartesian  coordinates  for  consistency. 
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Figure  14.  Surface  plots  of  a  Gaussian  beam  evolution  as  propagating  along  the  ;2- 
axis.  When  we  assume  z—lm,  (a)  is  an  intensity  profile  at  the  beam  waist  (2;=0.5m). 
At  (b)  (^=0.75m)  and  (c)  (^=lm),  we  can  see  the  expansion  of  the  beam  as  it  passes 
along  the  .s-axis  from  the  waist. 
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Let  us  assume  a  trial  solution  for  Equation  (IV. 7)  as 


V  =  S'  -  «  -  exp  ^  P[z)  + 
\wJ  \w/  [_  \ 


(IV.25) 


where  5-  is  a  function  of  x  and  2;  (because  w  depends  on  z),  and  h  is  a  function  of 
y  and  2;  to  account  for  diffraction.  Now,  we  should  find  g(x/w),  h(y/w),  P{z),  and 
q{z)  to  get  the  complete  solutions. 

For  convenience,  let  us  find  derivatives  of  g  and  h  with  respect  to  x  and  y.  By 
using  the  chain  rule,  we  get 


dg  dg  d{x/w)  1  dg  1  , 
dx  d{x/w)  dx  wd{x/w)  ' 

d‘^g  d  ( dg\  1  ( dg'\  1  dg'  d{x/w)  1 


dx"^  dx  \dxj  w  \  dx  J  wd{x/w)  dx 


(IV.26) 


1  \  dw  X  ,dtL 


dg  dg  d{x/w)  dg  f  1\  dw  x  ^dw 

dz  d{x/w)  dz  ^d{x/w)  \  w'^ )  dz  dz' 


(IV.27) 


where  (..)'  =  d{..)/d  [x/w).  We  can  easily  find  the  derivatives  of  h  through  similar 


processes: 


dh  1  I  d‘^h  1  n 

dy  w  ’  dy"^ 

dh  y  dw 

dz  w^  dz 


(IV.28) 


(IV.29) 


We  are  ready  to  take  derivatives  of  Equation  (IV.25).  After  some  tedious  work,  we 


obtain  the  following  results: 


^  hs"h  +  -!/h  ( -d  +  ah  (’Pi-ah('-x 


w  \q 


(IV.30) 


+  ^gh'  ( %)  +  sV*)  -  ah( 

dy^  w^  w  \  Q  /  \Q  /  \1 


(IV.31) 


=  e(-fr')  (-±^aA-±^ah'^y  +  aiJP-ah%^^)  .  (iv.32) 

dz  \  dz  w‘^  dz  dz  2q^  dz  J 
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Substituting  Equations  (IV. 30),  (IV. 31),  and  (IV. 32)  into  (IV.  7)  and  dividing  by 
gh  exp  f zP  +  ,  we  arrive  at  the  following  complicated  equation 


1  g"  2ikx  g' 


Wz  g'  1  h”  2iky  h' 


,  Wz  h' 

-r 
w‘^  n 


g  wq  g  g  w"^  h  wq  h 

' - v - '  ' - v - 

independent  oi  y  =  —Ui  independent  of  x  =  —U2 


[  — 


k^  k"^  dq 

q2  q2 


only  =  0 


2ik  dP 

q  oz 

" - V - ' 

no  r  =  Ml  +  M2 


(IV.33) 


In  order  for  this  equation  to  have  a  solution,  let  us  assume  that  the  ^-independent 
terms  are  —Ui  and  the  ^-independent  terms  are  —U2.  The  r^-dependent  term  should 
be  identically  zero  since  it  must  be  valid  for  an  any  location  in  the  transverse  plane 
[11].  In  the  end,  the  terms  left  at  the  end  should  be  ui  -\-U2  so  that  the  left-hand  side 
sums  to  zero. 

Let  us  consider  the  terms  of  Equation  (IV.33)  to  get 


,  k‘^  f  dq 


q^  \dz 


1=0. 


(IV.34) 


We  then  have  the  familiar  relation,  (dq/dz)  =  1,  which  we  have  seen  in  the  funda¬ 


mental  mode  derivation.  Thus,  we  have,  easily: 

1  _  1  ,  A 

q  P 


(IV.35) 


Next,  multiplying  the  a:-dependent  terms  of  Equation  (IV.33)  by  gw"^  and 
collecting  terms  in  the  derivatives  of  we  have 


//  ^ ,  f  "w  dw  \  ,  9  ^ 

g  +  2ikx - 9  +  uiw^g  =  0. 

\q  dz  J 


(IV.36) 


Let  us  calculate  only  the  2nd  term.  Recalling  Equations  (IV.  15)  and  (IV.  16)  and 
substituting  these  two  relations  into  the  2nd  term,  we  have 


,  f  w  iX  dw\  .  27r  /  iA 

2ikx  —  H - -^19  =  2.1X—  — 

\  K  'KW  dz  J  A  \  TTW 
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where  dwjdz  =  w/R.  We  can  easily  figure  out  this  relation  when  we  take  the 
derivative  of  w  with  respect  to  2;.  We  can  then  rewrite  Equation  (IV. 36)  as 


9 


// 


.  9 

— g  +uiw  ^  =  0, 
w 


(IV.37) 


which  is  close  to  a  form  of  the  Hermite  equation  y"  —  2xy'  +  2my  =  0  [16] .  In  order 
to  transform  Equation  (IV.37)  into  the  precise  form  of  a  Hermite  equation,  let 


X  = 


(IV.38) 


Then  we  have 


dg 

d  {x/w) 


Inserting  Equation  (IV. 39)  into  (IV.37)  yields 


d^g 


+  2mg 


0, 


(IV.39) 


(IV.40) 


which  is  exactly  the  same  as  the  Hermite  equation  where  2m  =  uiw‘^/2.  li  m  = 
0, 1,  2, . . . ,  then  a  solution  of  the  Hermite  equation  is  the  Hermite  polynomial  Hm{x) 
given  by  [16] 


six)  =  (X)  =  {bil)  =  (e-’^O  .  (IV.41) 

There  are  similar  solutions  for  the  y  dependent  terms  of  Equation  (IV. 33),  i.e.,  h{^)  — 
HniO  —  Hn  iy^y/w),  where  n  =  0, 1,  2, ... .  Some  of  the  Hermite  polynomials  of 
low-order  are  listed  in  Table  I. 

Now  we  are  left  with  only  the  last  terms  of  Equation  (IV. 33)  which  will  give 
a  solution  for  P{z).  According  to  Equation  (IV.40),  we  know  the  fact  that 


4m 


Ml 


U2 


4n 


By  using  these  two  relations,  we  rewrite  the  r- independent  terms  as 

2?A’  dP  4 

- 2k^  =ui  +  U2  =  ^{m  +  n).  (IV.42) 

q  oz 
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Ho{s)  =  1 
IIi(s)  =  2s 
H2{s)  =  4s2  -  2 
Hsis)  =  8s^  -  12s 
H4{s)  =  16s^  -  48s2  +  12 


Table  1.  Hermite  polynomials.  From  [16] 


Substituting  the  middle  terms  of  Equation  (IV.  18)  into  1/q,  we  obtain 


dP  i  2  .  ,  iz  ,  ,  N  ^0 

^  +  n)  =  — - r  -  (m  +  n  +  1)^- - r, 

oz  q  kw^  z^  +  Zq^  z^  +  Zq^ 

where  2/{kw‘^)  =  z^j^z^  +  Zq^).  After  integrating,  we  obtain 


i{m  +  n  +  1)  tan 


(IV.43) 


(IV.44) 


W  /  Zq 

which  is  similar  to  Equation  (IV. 22)  for  the  fundamental  mode.  Combining  all 
sources,  such  as  Equation  (IV. 41)  for  Equation  (IV.44)  for  iP  and  Equation  (IV. 35) 
for  1/g,  we  hnally  have  the  following  solution  for  higher-order  modes: 


/  - 
Ym.n 

W 


2-] 

Wy 


2-] 

w 


exp 


/y 

i\j !  /  V  _ 1/6 

—  (m  -|-  n  -|-  1)  tan  — 

ZjR  Zq 


Utilizing  this  solution,  Equation  (IV. 4)  can  be  written  in  the  form 

,  ^ 

2i'i  e-'’-'-" 


_^TT 

^m.n 

W 


2-]H„ 


w 


exp 


Amplitude 


(IV.45) 


\ 


kr^ 


kz  +  —  {m  +  n  +  1)  tan  — 

zR  Zq 

L  V  Phase  (j)  J 


(IV.46) 

where  the  positive  integers  m  and  n  identify  the  modes.  These  integer  indices  deter¬ 
mine  the  shape  of  the  beam  intensity  in  the  x  and  y  direction,  respectively.  Several 
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examples  of  higher-order  modes  are  shown  in  Figure  15.  The  intensity  maxima  are 
separated  by  anti- nodal  lines.  Note  that  mode  label  m  equals  the  number  of  intensity 
anti- nodes  proceeding  from  left  to  right,  and  n  the  number  proceeding  from  top  to 
bottom. 


- ►X  Wqj  M32 

Figure  15.  Intensity  patterns  for  Hermite  Gaussian  modes  in  the  transverse  plane,  x 
and  y. 
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V.  FEL  THEORY  WITH  HERMITE  GAUSSIAN 

MODES 


In  the  previous  chapter,  we  found  solutions  of  a  wave  equation  for  a  Hermite 
Gaussian  (HG)  optical  beam  in  free  space.  In  the  actual  FEL  operation,  a  beam  of 
electrons  interacts  with  the  Gaussian  optical  beam,  which  is  bounded  by  the  spherical 
resonator  mirrors,  as  shown  in  Figure  16. 


Figure  16.  A  schematic  of  the  FEL  oscillator  shows  the  Gaussian  beam  (blue),  with 
radius  w{z)  expanding  away  from  a  waist  wq  centered  in  the  undulator  (green)  of 
length  L.  Also,  the  optical  phase  (l){z)  (black)  changes  along  L.  One  electron  path, 
on-axis  at  r  =  0,  is  shown,  and  all  transverse  dimensions  are  exaggerated.  After  [24] 


This  chapter  explores  the  EEL  interaction,  considering  the  characteristics  of  HG  res¬ 
onator  modes.  We  will  hnd  a  modihed  pendulum  equation  describing  the  microscopic 
electron  evolution  and  determine  the  design  of  the  HG  mode  which  can  maximize  the 
interaction  with  the  electron  beam. 
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A.  MODIFIED  PENDULUM  EQUATION 

The  electric  field  at  the  mode  center  (r  =  0)  in  a  fundamental  Hermite  Gaus¬ 
sian  mode  as  derived  in  Chapter  IV  is 


£■(0,  z)  =  £o-^  exp  [i  {kz  +  (t){z))] , 
w[z) 


(V.l) 


where  Eq  is  the  electric  field  amplitude  at  the  mode  waist,  k  =  27:/ X  is  the  mode  wave 
number,  and  (/){z)  =  —  tan“^  [{z  —  L/2)/ Zq]  is  the  mode’s  phase.  Also  recall  that  the 
mode  radius  is  w{z)  =  Wq[1  -\-  {z  —  where  Zq  =  /X  is  the  Rayleigh 

length,  and  wq  is  the  mode  waist  radius  at  2;  =  L/2.  Although  we  will  only  deal 
with  the  fundamental  Gaussian  mode,  results  can  easily  be  extended  to  higher-order 
modes  [24],  The  mode  is  centered  at  =  1/2  (z  =  L/2),  since  this  point  is  near  a 
broad  maximum  in  gain  [24],  Using  dimensionless  time  r  =  z/L  =  ct/L  in  w{z)  and 
(/){z)  gives 


w‘^{t)  =  Wq^ 


(V.2) 


and 


^(r)  =  —  tan  ^ 


(V.3) 


Recall  the  pendulum  equation  derived  in  the  previous  chapter: 


00 

C  =  I'  =  a(T)  cos  [<  +  <I>(t)]  , 


(V.4) 


where  a(r)  is  the  dimensionless  field  amplitude  in  Equation  (V.l)  by  a^r)  =  aQWo/w{T), 
where  ao  =  AirNeKLEQ/'^'^mc^  is  the  dimensionless  optical  wave  amplitude.  Insert¬ 
ing  Equation  (V.3)  into  (V.4)  gives  us  the  modified  pendulum  equation,  describing 
an  electron’s  evolution  in  an  HG  mode.  The  parameter  L/ Zq  compares  the  undulator 
length  L  with  the  Rayleigh  length  Zq.  The  Gaussian  beam  radius  w  and  phase  (p 
depend  on  L/zq,  as  we  saw  in  Equation  (V.2)  and  (V.3).  In  order  to  illustrate  some 
of  the  effects  of  the  Gaussian  beam,  consider  L/ Zq  <^1  so  that  we  almost  have  plane 


waves  [24],  To  lowest  order  in  L/zq,  we  can  then  attain  the  following  relations  with 
L/zq  <  1: 

HJ  ^7” ^  ^  Wo,  (V.5) 


thus 


a(r)  =  ao 


Wq 

w{t) 


ao, 


(V.6) 


and 

using  the  small-angle  approximation.  Substituting  Equations  (V.6)  and  (V.7)  into 
(V.4)  yields  the  modified  pendulum  equation  for  the  case  of  low  gain  (j  <  tt) 


=1^ 


ao  cos 


C„  +  %r-|-(r-l 


ao  cos 
field 


new  phase  new  resonance 


(V.8) 


where  C  =  Co  +  t'o'r,  the  solution  of  C  to  zeroth-order  in  ao,  has  been  used  in  weak 
optical  fields  (uq  <  tt).  The  electrons  have  a  new  initial  phase.  But  the  shift  Co  ^  Co  + 
Lj^Zo  in  each  phase  is  not  significant  since  the  beam  is  uniformly  spread  over  each  op¬ 
tical  wavelength  [24].  The  gain  spectrum  has  the  form  [2  —  2  cos  vo*  —  t'o*(sini^o*)] 
derived  in  Chapter  III,  where  now  Vo*  =  Vo  —  L/ Zovasm  HG  mode.  The  gain  spectrum 
is  now  shifted  away  from  the  exact  resonance  by  an  amount  of  L/zo- 


B.  GAIN  IN  GAUSSIAN  MODES 

In  a  weak  optical  field  ao  <  tt  and  for  low  gain  j  ^  tt  ,  the  gain  in  the  HG 
modes  becomes 

G{t)  oc  jFa,  (V.9) 

where  j  is  the  dimensionless  current  density,  and  Fq  is  the  filling  factor  for  HG  modes. 
Now,  let  us  find  the  maximum  filling  factor  averaged  over  the  undulator  length  L  in 
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Consequently,  the  maximum  filling  factor,  ,  occurs  at  =  Ll\/\2  ~  {0.3)L. 

We  have  analytically  derived  that  the  gain  spectrum  is  shifted  by  pq  ^  pq  — 
L/ Zq  when  L/ Zq  <C  1.  We  have  also  shown  analytically  that  the  optimum  gain  occurs 
when  L/ Zq  ~  3.  Numerical  simulations  have  shown  that  the  shift  in  the  gain  spectrum 
still  occurs  for  values  oiLj zq  up  to  ~  10,  and  is  still  roughly  given  by  the  approximate 
shift  Pq—>Pq  —  L/zq  [26]. 
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VI.  HYBRID  (HG- WAVEGUIDE)  RESONATOR 

MODES 


For  short  wavelength  operation  of  an  FEL,  an  open  resonator  bounded  by 
spherical  reflecting  mirrors  is  feasible,  since  the  laser  beam  maintains  a  small  cross- 
sectional  area  as  it  evolves  along  the  interior  of  the  undulator.  However,  when  the 
FEL  operates  at  long  wavelengths,  such  as  in  the  sub-millimeter  or  THz  region,  the 
expansion  of  the  laser  beam  due  to  diffraction  can  be  a  problem,  as  shown  in  Eigure 
17. 


(d)  (e)  (f) 

Eigure  17.  Surface  plots  of  two  Gaussian  beams  propagating  along  the  ^-axis  without 
waveguides.  Compared  to  the  evolution  of  a  short  wavelength  beam  (Figures  (a),  (b) 
and  (c)  -  3/im),  the  evolution  of  a  long  wavelength  beam  (Figures  (d),  (e)  and  (f)  - 
30/im)  obviously  has  a  larger  amount  of  expansion. 


To  accommodate  lager  diffraction  at  longer  wavelengths,  we  need  to  have  a  larger 
undulator  gap  g,  but  the  larger  magnet  spacing  reduces  the  strength  of  the  magnetic 
field  according  io  B  =  1.95i?mac/ exp(— tt^t/Ao),  where  Aq  is  the  undulator  period  [2]. 
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Consequently,  the  gain  G  will  be  decreased  proportionally:  G  (x  B"^.  A  remedy  is  to 
use  a  waveguide  structure,  shown  in  Figure  18,  to  confine  that  the  radiation.  This 
chapter  analyzes  the  Hermite  Gaussian- Waveguide  resonator  modes,  which  will  be 
designated  as  Hybrid  resonator  modes.  The  modes  include  characteristics  of  both 
Gaussian  and  waveguide  modes.  In  order  to  understand  the  attributes  of  the  waveg¬ 
uide,  we  solve  the  wave  equation  and  explore  the  phase  and  group  velocities  of  the 
laser  pulse  propagating  through  the  waveguide. 


Figure  18.  Geometry  of  a  waveguide  resonator.  The  origin  is  located  at  the  center  of 
the  waveguide,  with  y  =  0  at  the  bottom  of  the  waveguide.  Curved  surfaces  represent 
the  mirrors.  From  [22] 


A.  SOLUTION  OF  WAVE  EQUATION 

The  Hybrid  resonator  modes  have  been  studied  in  detail  by  Elias  [17],  Degnan 
[18],  and  others.  This  section  follows  Elias’s  analysis,  which  is  the  basis  of  the  UCSB 
FEE.  To  find  solutions  of  a  wave  equation,  assume  that  the  electrical  conductivity  of 
the  walls  in  Eigure  18  are  infinite.  Eirst  find  solutions  to  the  wave  equation  inside  a 
perfect,  infinitely- long,  rectangular  waveguide  in  terms  of  Gaussian  waves  propagating 
along  the  z-axis.  These  manipulations  are  similar  to  the  derivations  of  Chapter  IV; 
the  difference  is  that  we  need  to  consider  waveguide  properties  in  the  r/-direction. 
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1.  Fundamental  Mode 

Each  cartesian  component  of  the  electric  (or  magnetic)  field  must  be  satisfied 
by  the  Helmholz  equation  presented  in  Chapter  IV: 


{V^  +  k^)u{x,y,z)  =  0,  (VI.l) 

where  k  =  ujjc  =  27r/A,  and  a  harmonic  time  variation,  has  been  assumed. 

Note  that  u{x,  y,  z)  can  represent  either  a  component  of  the  electromagnetic  field  or 
of  the  vector  potential  of  the  light.  We  are  interested  in  finding  solutions  describing 
Gaussian  waves  propagating  along  the  .s-axis  without  phase  variations  along  the  y-axis 
[17].  We  express  u{x,y,z),  an  electromagnetic  wave  traveling  along  the  2;-direction 
within  the  waveguide,  as 

u{x,  y,  z)  =  '4){x,  sin  (mry/b) ,  (VI. 2) 


where  '0(x,  z)  is  a  slowly-varying  complex  function,  is  a  fast-varying  carrier 

wave,  sm{m:y /b),  b  is  the  distance  between  the  top  and  bottom  plates,  and  n  — 
±1,  ±2,  ±3, ...  is  an  integer  in  order  to  match  boundary  conditions  at  the  walls  of 
the  waveguide  (note  y  =  0  at  the  bottom  wall).  Inserting  Equation  (VI. 2)  into  the 
scalar  wave  equation  and  dividing  by  sin(n7ry/6)  yields 


d'^ip  d^ip 


dx"^ 


2ik  ^ 

OZ^  OZ 


=  0. 


From  this  equation,  we  find  that 


k„  =  \k^- 


ri^TT^ 

1^' 


(VL3) 


(VI.4) 


When  we  apply  the  paraxial  approximation  [12],  d'^'i/^ldz'^  <C  d‘^ip/dx‘^,  and  we  arrive 
at  the  paraxial  wave  equation 


8“^ 'll) 


2iki 


dip 


dx"^  '  dz 


Let  us  assume  a  trial  solution  of  the  form 


=  0. 


ip{x,  z)  =  exp 


P(z) 


2q(z) 


X 


(VI.5) 


(VI.6) 
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similar  to  the  Hermite-Gaussian  case,  except  that  here  there  is  no  y  variable.  We  need 
to  find  the  unknown  functions  P{z)  and  q(z)  in  a  similar  method  to  that  described 
in  Chapter  IV.  After  inserting  Equation  (VI. 6)  into  (VI. 5)  and  comparing  terms  of 
equal  power  in  x,  we  obtain  the  following  relations 


dq 

dz 


=  1  , 


(VI.7) 


(VI.8) 


and 

dP  _  i 
dz  2q  ’ 

In  order  to  find  P{z)  and  q{z),  we  use  the  definitions  of  w{z)  and  R{z)  again,  but 
these  are  not  quite  the  same  as  the  previous  parameters  since  a  new  Rayleigh  length, 
zr  =  WQ^kn/2  [17],  differs  from  the  previous  one  in  Chapter  IV,  Zq  =  rvoq^/X  = 
wo^k/2.  Using  this  new  definition  of  zr,  we  integrate  Equation  (VI.7),  yielding 


q  =  z-\-qo  =  z-  i 


.  ^0 


=  Z-  IZr, 


(VI.9) 


where  the  complex  beam  parameter  at  the  waist  is  qo  —  —iwo^kn/2.  Starting  with 
Equation  (VI.9),  q  =  z—izR,  and  employing  Equations  (IV. 15)  and  (IV. 16)  in  Chapter 
IV,  we  obtain  the  following  relation 


2q 


2R  ^  m;2- 


(VI.  10) 


To  determine  the  P{z),  we  insert  Equation  (VI.9)  into  (VI.8)  and  integrate  the  result, 
and  we  then  have 


iP{z)  —  In 


w 


Wq 


n 

2 


-tan 


-1 


Zr 


(VI.ll) 


using  the  definition  of  w{z).  Substituting  Equations  (VI. 10)  and  (VI.ll)  into  the  trial 
solution,  Equation  (VI. 6),  we  can  finally  get  the  solution  of  the  fundamental  mode 
inside  the  waveguide 


i){x,z) 


(VI.12) 
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This  solution  for  the  fundamental  mode  is  also  similar  to  the  Equation  (IV. 24)  in 
Chapter  IV,  but  now  kn  appears  instead  of  k,  and  the  Gaussian  form  is  only  in  the 
x-direction.  Using  this  solution,  Equation  (VL2)  can  be  written  in  the  form 


(  \ 

/  \ 

u[x^y^z)  =  \  — e  '  exp 

V  w 

% 

,  knx'^  1  1  2; 

knZ  H - ;; - tan  — 

2R  2  zr 

Amplitude 

Sr 

V  Phase  (p  / 

'  riTT 

sin(— y).  (VI.  13) 


2.  Higher-Order  Modes 

The  higher-order  solutions  in  the  waveguide  also  take  the  form  of  Hermite- 
Gaussian  functions  in  Gartesian  coordinates,  as  we  now  show.  Assume  a  trial  solution 
to  Equation  (VI. 5)  in  the  form 


'ip  =  g  exp 


T,  ,  2 

P+—X 

2q 


(VI.  14) 


where  is  a  function  of  x/w.  We  can  easily  find  the  derivatives  of  the  unknown 
function  g  when  we  follow  the  notation  mentioned  in  Chapter  IV.  They  are 


dx^ 


1  //  2  f ikji  \  f  f f kfi 

—g  +  -  — X  U+  —  ]  g  -  [—X  g 

w  \  q  j  \  q  J  \q 


w 


and 


dz 


1  idw  ,dP  iknx‘^  dq' 
-:^g  +  gi—  -  9^^-^ J  , 


(VI.15) 


(VI.  16) 


where  (..)'  =  d{..)/d{x/w).  Inserting  Equations  (VI.15)  and  (VI. 16)  into  (VI. 5)  and 
dividing  by  5" exp  (iP  +  iknx‘^/2q),  we  arrive  at  the  following  result 


1  g"  2iknxg'  dw/dzg'  ^ 

r)  ~r  ZtrvjiX  ^  X 

g  wq  g  g 

^  ^  ^ 


kr] 


kn^  dq 


q^  dz 


ikn 


dP 


dz 


X  only  =  ui 


x^  only  =  0 


no  X  =  ui 


(VI.  17) 

Utilizing  the  same  notation  as  in  Chapter  IV  allows  us  to  have  g{x/w),  P{z),  and 
q(z)  such  that 

g(-)=H^(-),  m  =  0,l,2,...,  (VI.18) 

\wy  \w/ 


45 


iP  =  In 


(VI.  19) 


in 

w  ^  2 

Wo 


—  i  {  m  —  1  tan  ^  — , 
2  /  zo 


and 


1  1  2i 

-  =  ^  + 


r.  •  ,  o  •  (VI.20) 

q  R  knW^ 

Substituting  these  three  equations  into  Equation  (VI.  14),  we  finally  have  the  following 
solution  for  higher-order  modes 


iWo 


'4^m,n  \l  Hffi  I  )  e 
w  \w 


X 


exp 


knX'^  (  1\  _i  .2 


(VI.21) 


Utilizing  this  solution,  we  can  express  the  optical  beam.  Equation  (VI. 2),  as 

/  \ 


w  \w 


exp 


Amplitude 


h 

k^z  - 

”  2i? 


m  +  -  1  tan 


-1 


zr 


\ 


Phase 


.  /"mr 
sm  \^—y 

(VI.22) 
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Figure  19.  Intensity  patterns  E'^rnn  for  Hybrid  modes  in  the  transverse  plane,  x  and 

y- 
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Note  that  n,  the  waveguide  constant  in  the  sine  function,  differs  from  the  previous 
n  of  Hn  in  Chapter  IV.  However,  this  n  together  with  m  determines  the  shape  of 
the  mode,  as  shown  in  Figure  19.  We  chose  the  sm(mry/b)  as  a  waveguide  term  in 
Equation  (VI. 2)  to  satisfy  the  boundary  conditions;  there  are  no  helds  for  n  =  0. 

3.  Solutions  for  E’^^nn  Modes 

Elias  and  Gallardo  [17]  proposed  two  sets  of  linearly  independent  solutions 
differing  by  the  choice  of  sm{mry/b)  or  cos{n7ry/b).  In  the  solutions  for  sm{n7ry/b), 
which  will  be  designated  as  modes,  the  dominant  electric  fields  are  parallel  to 

the  wide  dimension  of  the  waveguide,  as  shown  in  Figure  20. 


Figure  20.  Electric  field  lines  at  the  waist  of  the  waveguide,  which  shows  the 
mode. 

The  solutions  for  the  choice  cos{mry /b),  which  will  be  designated  as  modes, 

the  dominant  electric  fields  are  those  which  are  perpendicular  to  the  x-direction  of 
the  waveguide  not  shown  in  Eigure  20.  When  we  chose  the  E'^rnn  modes,  waveguide 
losses  are  negligible  since  the  electric  fields  near  the  plates  are  close  to  zero.  In  an 
FEE,  the  undulator  magnets  would  be  placed  along  the  y-sods,  so  that  an  electron 
wiggles  along  the  x-axis  and  couples  to  the  electric  fields  of  the  mode,  as  shown  in 
Figure  21. 
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Waveguide 


Figure  21.  Schematic  of  a  Hybrid  resonator,  showing  the  location  of  the  critical  com¬ 
ponents  after  [19].  The  hybrid  waveguide  consists  of  a  set  of  two  parallel  conducting 
plates  separated  by  a  distance  b.  A  relativistic  electron  wiggles  along  the  x-axis.  A 
hole  may  be  placed  in  the  middle  of  the  right  mirror  to  allow  output  of  the  laser  beam 
(KAERI  is  using  this  type  of  mirror  for  their  FIR  FEL  [20]).  Otherwise,  a  partially 
transparent  mirror  will  be  used  for  the  out-coupling.  In  case  of  entry  and  exit  of 
the  electron  beam,  a  chicane  is  usually  used  to  bend  the  electron  beam  around  the 
mirrors. 


We  are  now  ready  to  explore  the  solutions  for  modes.  Recall  the  pro¬ 

posed  solution  u(x,  y,  z)  =  ip(x,  z)  sin(n7r|//6)e*^"^.  We  can  then  express  the  dominant 
electric  fields  as 

(HIT  \ 

—yj  (VI. 23) 

where  Aq  is  a  constant.  Using  this  equation  and  Maxwell’s  equation  V  ■  E  =  0  in  free 
space,  we  can  obtain  the  Ez  component: 


E, 


Aq  dlpznn 
ikn  dx 


sin 


giknZ 


(VI.24) 


Note  that  the  component  is  much  greater  than  Ey  and  Ez,  i.e.,  E^  ^  Ey,  Ez-  Also, 
=  0,  as  represented  by  Eigure  20.  In  order  to  determine  the  magnetic  fields,  we 
use  Maxwell’s  equation  V  x  E  =  —dB/cdt.  When  we  consider  only  the  y  component 
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of  this  equation,  we  have 


cdt 


dE, 

dz 


-{ikn)Ao^rnn  Sm 


JknZ 


and 


By  ^  ^Ao^mn  Sin  (^yj  (VI.25) 

where  d^l^mn/dz  has  been  ignored,  since  the  amplitude  of  the  longitudinal  components 
of  the  fields  are  much  smaller  than  the  corresponding  transverse  components,  and  the 
harmonic  variation  has  been  assumed,  respectively.  Next,  let  us  discuss  the  .2 
component  of  V  x  E  =  —dG/cdt.  We  then  have 


dB,  dE^ 


cdt 


dy 


\  ik  z 

—  Ao'lprnn  COS  — yj  6  , 


and 

Bz  ~  cos  (^y)  (VI.26) 

where  a  harmonic  variation  has  been  assumed  again.  Lastly,  the  boundary 

conditions  applied  on  the  E’^^nn  modes  are  E^^E^  =  0  at  y  =  0  and  y  =  b.  A 
summary  of  the  E'^rnn  modes  solutions  is  given  by  the  following  table. 


E^  Ao'ipmnAn{mvy/b)e^^’^^ 

Ey  =  0 

Ez^  -  {AQ/ikn)  {diprnnidx)  sin(n7ry/6)e*^’»"^ 


Table  2.  E'^rnn  -  modes  solutions 


Starting  with  these  solutions,  we  hnd  a  relation  between  the  coefficient  Aq  and 
the  power  of  the  electromagnetic  held.  The  expression  for  total  energy  density  [21] 
is: 

u=^{E^  +  B^).  (VI.27) 
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B.  PHASE  AND  GROUP  VELOCITY 

The  general  expressions  for  phase  and  group  velocities  of  a  electromagnetic 
wave  are  Vp  =  u/k,  and  Vg  =  du)/  dfc.Then  from  Equation  (VI. 4)  with  <C  zr,  we 
have: 

kn^  1  1  /n7r\2 

Vpn'^  (?  Vcnfe/  ' 

Therefore, 

Up,  =  ""  (VI.32) 

Y  1  —  {n'Kcluh'f 

Note  that  the  phase  velocity  of  the  wave  depends  on  the  waveguide  number,  n  = 
1,2,3,...,  the  frequency,  cu,  and  the  height  of  the  waveguide,  h.  Also,  the  phase 
velocity,  Up,,  is  faster  than  the  speed  of  light,  c,  since 
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The  group  velocity  of  the  wave  can  be  derived  from  Equation  (VI. 4).  Differentiating 
this  equation  with  respect  to  cn,  we  get 


2kn 


d/^7x 

den 


u 


=  2— 


dkr, 


u 


den  knC^ 


Hence,  the  group  velocity  is 


den  k'ic  (? 


Vgn  — 


dikr, 


en  V, 


pn 


We  can  rewrite  Equations  (VI. 43)  and  (VI. 45)  as 
c  c 


"^pn 


1  —  (mrc/ijby 


1  —  {cnx/2xfhy 


1  -  {n\l2hy 


and 


^9n-  -C'\/l  (25)- 


~^pn 


Thns,  VpnVgn  =  (?  is  also  satisfied  inside  the  wavegnide. 


(VI.33) 


(VI.34) 


(VI.35) 


(VI.36) 


Note  that  when  the  wavelength  is  equal  to  the  cut-off  wavelength,  i.e.,  Ac  =  26 
and  9  —  90°  as  shown  in  Eigure  22,  the  phase  velocity  is  infinite  and  the  group  velocity 
is  zero  so  that  no  energy  propagates  along  the  wavegnide. 


Figure  22.  A  sinusoidal  electromagnetic  wave  entering  one  end  of  the  wavegnide  with 
its  electric  field  parallel  to  the  wide  dimension  and  propagating  at  an  angle  9  will  be 
successively  reflected  by  the  top  and  bottom  walls  and  follow  a  zig-zag  path.  A  is 
near  the  cut-off  wavelength.  From  [23] 
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The  cut-off  frequency  fc  is  given  by 


52 


VII. 


FEL  THEORY  WITH  HYBRID 
(HG- WAVEGUIDE)  MODES 


A  Hybrid  mode  occurs  in  a  waveguide  consisting  of  two  flat  conducting  plates 
seperated  by  a  distance  6,  each  parallel  to  the  x-z  plane.  This  waveguide  serves  as 
an  optical  beam  pipe  between  the  resonator  mirrors.  In  principle,  there  are  image 
charges  in  the  conducting  walls  that  may  affect  the  electron  trajectories,  but  they 
are  very  small.  More  important  things  are  the  space  charge  forces  in  the  beam  itself 
which  are  assumed  negligible  here,  as  in  most  FELs.  In  this  chapter,  we  will  explore 
a  modifled  pendulum  equation  appropriate  to  Hybrid  modes.  We  will  then  design  the 
relevant  FEL  which  results  in  maximum  energy  gain,  using  the  same  analysis  as  in 
Chapter  V. 


A.  MODIFIED  PENDULUM  EQUATION 


The  electric  field  at  the  mode  center  (a:  =  0,  r/  =  6/2)  in  a  fundamental  Hybrid 
mode  as  derived  in  Chapter  VI  is 

^  ^o^^^exp  [i  {kiz  +  (j){z))] ,  (VH.l) 

where  Eq  is  the  electric  held  amplitude  at  the  mode  waist,  ki  =  ^/k"^  —  and 

(p{z)  —  —Q.bicai~^[{z  —  L/2)/zB\,  where  zr  =  wo^ki/2  is  the  Rayleigh  length  for 
this  Hybrid  mode  case,  and  Wq  is  the  mode  waist  radius  at  .2  =  L/2.  Although  we 
will  only  deal  with  the  fundamental  Hybrid  mode,  results  can  easily  be  extended  to 
higher-order  modes.  We  can  again  introduce  the  dimensionless  time,  r  =  z/L,  so 
that  the  mode  is  centered  at  =  1/2,  near  a  broad  maximum  in  gain  [24].  Through 
a  derivation  similar  to  that  in  Chapter  V,  we  have 

and 


(/)(r) 


(VIL3) 
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Our  original  pendulum  equation  in  Chapter  III  is 

oo 

(  =  u=  a(T)  cos  [C  +  (pir)] ,  (VIL4) 

where  a(r)  is  the  amplitude  factor  in  Equation  (VII.  1),  where  ao  =  Atv N eK LEq / 
is  the  dimensionless  optical  wave  amplitude.  Using  the  approximation,  L/zr  1  de¬ 
scribing  a  near  plane  wave,  gives  us  to  lowest  order  in  LjzR 

w{t)  ~  tco,  (VII. 5) 


so  that 


a(r)  =  ao 


ao, 


(VII.6) 


and 

Substituting  Equation  (VII.6)  and  (VII. 7)  into  (VII. 4)  yields  the  modified  pendulum 
equation  for  the  case  of  low  gain  (j  <  tt),  derived  here  for  the  first  time: 


^ =1/  ao  cos 


ao  cos 
field 


Co 


^Zr 


+  T  \  Vo  - 


new  phase 


‘^Zr 

new  resonance 


(VII.8) 


where  C  =  C,o  +  vot  has  been  used  assuming  weak  optical  fields  (ao  <  tt).  The  electron 
has  a  new  initial  phase  (poE L/2zr^  but  this  shift  can  be  ignored  because  the  electrons 
are  uniformly  spread  along  each  optical  wavelength.  However,  we  cannot  ignore  the 
new  resonance  parameter,  vq  (z/q  —  L/2zr),  which  causes  a  measurable  change  in 
the  EEL  operation.  In  a  weak  optical  field  ao  <  vr,  the  gain  spectrum  has  the  form 
[2  —  2  cos  1^0*  —  z^o*(sinz/o*)]  /vq*^  derived  in  Chapter  III,  where  now  i^o*  =  vq  —  L/2zr 
in  the  Hybrid  mode. 
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B.  GAIN  IN  HYBRID  MODES 


In  a  fundamental  Hybrid  mode,  the  transverse  mode  area  can  be  used  to 
calculate  the  filling  factor.  This  affects  the  gain,  since  gain  is  proportional  to  the 
filling  factor  {G  oc  JFh).  The  cross-section  area  of  the  elliptical  Hybrid  mode  can  be 
expressed  irr^Vy,  where  =  wi^z),  Vy  ~  6/2,  and  6  is  the  height  of  the  waveguide  as 
shown  in  Figure  23. 


Figure  23.  Cross-section  of  the  fundamental  Hybrid  mode.  The  area  of  the  elliptical 
cross  section  of  the  mode  is  TTr^Ty  ~  7rw{z)b/2. 


By  using  a  similar  analysis  as  in  Chapter  V,  we  can  find  what  zr  value  gives 
the  maximum  Fr,  which  leads  to  the  maximum  gain.  In  a  weak  optical  field  ao  <  vr, 
the  gain  for  the  fundamental  Hybrid  mode  is  given  by 

G{t)  oc  JFh,  (VH.9) 


where  j  is  the  dimensionless  current  density,  and  Fr  is  the  new  filling  factor  for  the 
Hybrid  mode,  which  we  will  explore  in  detail  through  this  section.  The  averaged 
filling  factor  Fr  becomes 


electron  beam  area 


2'Krh^ 


^  (VII.  10) 

optical  beam  Area  7rbw{z) 

where  is  the  electron  beam  radius  (assumed  constant),  w{z)  is  the  optical  mode 

semi- major  axis  along  x,  and  6/2  is  the  optical  mode  semi- minor  axis  along  y.  We 

can  write  the  area  of  the  fundamental  Hybrid  mode  as 


.  ,  ,  irb  ,  .  7r6  /  ^ 

A{z)  =  —w{z)  =  —Wo  1  -I - ^ 

Zr^ 


1/2 


(VH.ll) 
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where  zr  =  Wo^ki/2,  and  ki  =  {k‘^  —  Thus,  the  averaged  mode  area  be¬ 

comes 


TrbwoZR 

2L 


7rb 

y/WiL 


A{z)  dz  = 


2zr  \4:Zr‘^ 

2  U;2ii 


\ 

j  +  zr^^^  In 


(VII.  12) 


where  wq  =  {2zR/kiy^‘^  has  been  used  in  the  last  step.  We  need  to  find  the  zr'^^^^ 
value  which  leads  to  the  maximum  filling  factor,  and  minimum  A.  However,  dA/  dzR 
is  a  complicated,  so  let  us  find  zr^^^  numerically.  From  the  numerical  solution,  we 
obtain  the  estimated  zj^^  value,  which  yields  the  smallest  optical  area,  as  shown  in 
Figure  24. 


Figure  24.  Averaged  optical  mode  area  A  for  fundamental  Hybrid  mode  versus 
Rayleigh  length  zr,  when  b  =  1cm  and  L  =  Im.  We  see  from  this  graph  that 
occurs  at  zr^^'^  ~  (0.4)L. 
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Therefore,  the  maximum  filling  factor,  occurs  at  ~  (0.4)L,  giving  a  peak 

value  of  G/ jFn  (gain  normalized  to  j  and  Fh)  of  about  0.135  and  at  ~  3.85. 
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VIII. 


CONCLUSION 


In  this  thesis,  we  have  reviewed  the  basic  PEL  theory  and  found  solutions  of 
the  wave  equation  for  Hermite  Gaussian  (HG)  resonator  modes  in  free  space.  We 
have  then  analyzed  the  effects  of  HG  modes  on  EEL  theory,  including  the  pendulum 
equation  and  gain.  We  have  then  studied  the  Hybrid  (HG- Waveguide)  resonator 
modes.  There  are  crucial  for  long  EEL  wavelength,  where  a  waveguide  must  be  used 
in  the  resonator  cavity  to  prevent  the  optical  beam  from  diffracting  into  contact  with 
the  undulator.  We  have  analyzed,  for  the  first  time,  EEL  operation  with  these  Hybrid 
modes. 

A  waveguide  in  the  resonator  cavity  produces  a  change  in  the  EEL  resonance 
parameter.  This  change  has  the  effect  of  shifting  the  gain  spectrum  away  from  the 
traditional  (free  space)  value  of  resonance.  In  a  fundamental  Hybrid  mode,  the  trans¬ 
verse  mode  area  must  be  used  to  compute  the  filling  factor,  which  also  affects  the 
gain.  We  have  found  a  new  optimized  Rayleigh  length,  associated  with  the  waveguide, 
which  allows  the  EEL  to  have  maximum  gain.  Einally,  we  recommend  the  following 
future  work: 

•  Examine  the  beam  coupling  to  free  space  mode  outside  the  Hybrid  mode 
region. 

•  Examine  the  output  coupling  with  a  partially  transmissive  mirror  versus  beam 
output  that  using  a  mirror  with  a  hole. 
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